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1 Introduction: Fluid models of nuclei 


Fluid models of nuclei have been used for quite a long time. The time-honored 
liquid drop model Jl|| is extremely successful in describing the ground state 
bulk properties of nuclei by utilizing extremely simple assumptions. These 
assumptions are based on the observation of the saturation phenomenon in 
nuclei as the predominant effect in determining the binding energies. The 
analogy of the nucleus to a fluid is also borne out in the almost constant 
density of nuclei throughout the periodic table and the lack of localization of 
nucleons around fixed positions, as it would be for a solid, due to zero point 
motion. The nucleus is a quantum fluid and cannot be modeled as a usual 
fluid for which the molecules behave to a large extent classically. 

More recently, the need to include relativistic effects in the description of 
hadron dynamics inside the nucleus lead to the development of the nuclear 
mean field models |2j]. These theories can be viewed as fluid models in which 
the mean field plays the role of the continuum (hydrodynamical regime), 
while the elementary fermion fields are their microscopic source. Again, 
the model has proven fruitful in the prediction of ground states as well as 
collective states of nuclei. 

On the other hand, nuclei -especially light ones- possess shell structures 
that are usually traced back to the single particle energy levels. This is 
clearly seen in the oscillations of nuclear binding energies as a function of 
mass, as well as in the existence or inexistence of stable nuclides. 

Heavy nuclei, however, seem to behave as deformed rigid bodies. This 
collective behavior is put in evidence by the large quadrupole excitation 
amplitudes. 

The question of the nature of the state of matter in the nucleus conse¬ 
quently depends on both the excitation energy, the mass and the density of 
the nucleus. 

In any event, the nucleus cannot be compared to a solid. Even heavy 
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nuclei have a rather soft incompressibility factor. Being compressible, they 
resemble more a gas or a compressible liquid, as advocated by the propug- 
nators of the liquid drop model. 

If one agrees with this view of a nucleus as a fluid, liquid or gas depending 
on the temperature, then it is valid to consider its description in terms of a 
few degrees of freedom as it is done for conventional fluids. The quantum 
nature of the nucleons demands the treatment to be quantal. The high speeds 
attainable by nucleons in the nucleus requires a relativistic treatment. 

One key ingredient that is absent in all the hydrodynamical models is the 
nucleon finite size. The charge radius of the proton is not negligible com¬ 
pared to nuclear dimensions, in contradistinction to molecules in a liquid. In 
the nonrelativistic Hartree-Fock or relativistic mean field Hartree treatments 
of nuclei, the nucleon is treated as pointlike. Some information about the 
finite extent of the nucleon is sometimes included in the form of short range 
correlations, short range hard-core, etc. 

The purpose of the present work is to try to deal in a coherent manner 
with the fluid aspect of nuclei and the finite size of the nucleon. 

Much has been learned in the past 40 years about the nucleon. Despite the 
fact that QCD is the accepted theory of strong interactions, it is practically 
impossible to use it as a tool to generate baryons as confined objects of quarks 
and gluons, and even less to describe nuclei. We here resort to a low energy 
model of strong interactions, namely, the Skyrrne model ||. The issue of finite 
nucleon size is dealt with by describing the nucleon as a topological soliton, 
the skyrmion. The fluid aspect is addressed by using a relativistic mean 
field theory in the spirit of the Walecka model ||. The scheme is simplified 
by resorting to a dilute fluid approximation to be explained below. In this 
manner the nucleons are essentially free, interactions being mediated by the 
mean fields. Section 2 summarizes the formalism and previous results both 
for nuclear matter and finite nuclei. Section 3 is devoted to the energy levels 
of a skyrmion in the nucleus. Section 4 treats the effects of the meson fields 
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on the properties of the skyrmion. 


2 A skyrmion fluid model of nuclei 

In the early 60’s Skyrme developed a topological soliton model for nucleons^]. 
In this model, baryons emerge as classical topological solutions of a nonlinear 
meson lagrangian. More than a decade ago this perspective was taken up 
by Wittenp. He showed that baryons may indeed appear as solitons in the 
large N c limit of QCD. The development of collective coordinate quantization 
of soliton rotations || gave a tremendous boost to the held. The door was 
now open for the investigation of baryons and nuclei. The Skyrme model and 
its topological solitons: skyrmions, may be regarded as the effective degrees 
of freedom of low energy baryon physics. This is especially appealing due to 
the difficulty of solving QCD explicitly. 

The Skyrme model has had moderate success in dealing with the nucleon- 
nucleon interaction, including the attractive isoscalar central potential for 
which several mechanisms were proposed |5J. It is difficult to solve the model 
exactly for large baryon number structures. Static solutions of various ge¬ 
ometries for some special baryon number cases were recently found using the 
rational map ansatzjTj. These solutions are of a limited practical use for 
the treatment of nuclei. Nuclei appear to be extremely dynamical assembly 
of nucleons. The special static solutions found in the literature resemble a 
crystal P|. Some of these solutions are in conflict with the spin-isospin as¬ 
signment of nuclear ground states[|j]. We here opt for an approach in which 
the nucleus is viewed as a dynamical aggregate of nucleons. For this purpose 
it is convenient to follow the path of fluid-like mean held theories, such as 
the relativistic model of Waleckap]. 

In the relativistic mean held theory of Walecka the nucleons are de¬ 
scribed using Dirac wave functions. The nucleons interact with meson helds. 
The meson helds are in turn determined by the various densities of the 
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baryons(scalar, vector,etc.). The nucleons have no internal structure. They 
are considered elementary pointlike objects. The model is simple enough 
to handle the basic features of nuclei and nuclear matter, and yet powerful 
enough to predict ground states and excited states with a relatively small 
amount of parameters that enter in the meson self-interactions and meson- 
nucleon coupling constants. The model is used at the first quantized level. 
Quantum corrections spoil the predictions. Being only a low energy effective 
model for nuclei, it is valid to consider it within its realm of applicability: 
the mean field level. 

The same path will be followed here. We will construct and use a mean 
field model of nucleons in terms of skyrmions, instead of Dirac pointlike ob¬ 
jects. The task becomes more complicated due to the finite volume occupied 
by the skyrmion. In some sense the Dirac wave function represents only 
the center of mass dynamics, while, the first quantized skyrmion will carry 
information about the interior of the nucleon. 

The skyrmion lagrangian is built by demanding both isospin and chiral 
symmetries (SU(2)xSU(2)). The many-body model developed here still obeys 
the same principles. 

The mean fields considered in the present work are the minimal set in 
order to be able to describe low energy nuclear phenomenology. A scalar 
held provides the attractive interaction and binds nuclei. The coupling of 
the scalar to the soliton is taken from phenomelogically succesful models [ITdj 

tween skyrmions ||13||, that is absent for potentials generated from the product 
ansatz. 

In order to stabilize the nucleus and reproduce the saturation property 
we also need a repulsive interaction. We here follow the Walecka models and 
use a meson held that is an isoscalar-vector held, the u; meson. 

Consider a held theory lagrangian of skyrmions, a scalar ifj and the u 


11 


The added scalar brings about the intermediate range attraction be- 
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meson H 


C' — -^-'2 t/j T £-2 T £-& ^interaction D(0) T £-uj 

= 7^0 < 9 M 0 - g i , l --^tv(L^L fl ) + y^tr[£„, L0 2 - g v u;^ 

- + (!) 

Here 

L m = (2) 

where U(r,t) is the chiral held, F n is the pion decay constant and e the 
Skyrme parameter, 0 is the scalar, 7 , an energy scale-parameter, and, g t i ,, 
a coupling constant. We have omitted the p meson as it is our intention to 
treat symmetric nuclei only. 

In the present work we will not specify the potential of the scalar. It must 
contain a mass-term and higher order terms. 

Due to the fact that the sign of 0 is irrelevant, it leads to an attractive 
interaction for both positive and negative signs, we can redefine the held in 
the form 0 = 7 e a . We will henceforth absorb the coupling constant g^ in 
the definition of the constant F^, that is in any event a htted parameter in 
the Skyrme model. 

Similarly to the mean held theory of pointlike baryons, we consider an 
ensemble of essentially free skyrmions. Each skyrmion will be accompanied 
by its own scalar held and u helds. Although we deal with free skyrmions, the 
average properties of the ensemble are still included in distribution functions 
that depend on the density and temperature. 

The skyrmion is a topological soliton. The winding number of the soliton 
was already identihed by Skyrme with the baryon charge. It is absolutely 
conserved regardless of the dynamics. It is a geometrical property, not a 
Noether charge. We have to take care of this aspect when building many 
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body ansatze for the nucleon fluid. The baryon number should be exactly 
conserved. 

An ansatz that conserves baryon number exactly and still allows for a 
reasonable treatment of the dynamics of N skyrmions is the product ansatz 

Ub=n{?i Ri, Rn ■ ■ • j R-at) = U( r — Ri) U( r — R 2 ) • ■ • U(r — Rat), (3) 

Where Ri is the location of the center of skyrmion i. 

This ansatz is not an exact solution of the skyrmion sector of the problem. 
We mentioned static exact solutions 0 as being appropriate for a crystal-like 
structure of nuclei. It is our intention to treat the dynamics of the nucleons 
in the nucleus in a reasonable manner without sacrificing nucleon motion, 
hence we opt for this crude approximation. Despite being a rough ansatz, 
it is very maleable and quite appropriate for the nuclear case in which the 
nucleons are not tightly packed. Nucleons in the nucleus are separated by an 
average distance in the order of 2 fin. The solitons overlap at the ’surface’ 
of each skyrmion. This interaction is replaced by the effective meson fields. 
Instead of dealing with the complicated situation of the N-body interaction 
in all its aspects, a small set of mesonic degrees of freedom is chosen and the 
skyrmions interact solely with them. It is expected that the introduction of 
the mesons in the mean fields will compensate, at a phenomenological level, 
the lack of accuracy induced by the rudimentary product ansatz. 

For the meson fields we use (see section 4 for an improvement upon this 
approximation) 


04 T 0"2 T • • • T <tn, 


(7$ 5<j\ -\- 5@2 

(4) 

Oil + U>2 + • • • + UJn, 


Uq + <5 011 + S 0 J 2 + • • • + 

(5) 
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Where <Jo,uJo are the mean held constant values of the scalar and the 
u and da, du> represent the huctuations. In thermal equilibrium, the mean 
held helds will depend on the temperature T and the chemical potential /j. 
Self-consistency then demands the values of a 0 and uj 0 to be determined by 
the properties of the ensemble. In equilibrium the huctuations will be taken 
to vanish in accordance with the mean held approximation. 

The topological baryon density 


B^ = 




-tr 


247T 2 

with the the product ansatz of eq. 


( U ] d a U ) ( U'dpU ) (rfd^U 


gives 


( 6 ) 


B 0 — b\ + 62 + • • • + bjy, 

(7) 

The dynamics of the nucleons is that of a rigid soliton. No internal 
excitations are considered. Also skyrmion distortions will not be allowed, 
except for swelling and shrinking of the soliton in the medium. Only the bulk 
parameters of the skyrmion, mass and moment of inertia will be allowed to 
vary due to the presence of the mesonic medium. 

This is the basis for the dilute fluid approximation. It is related to the 
mean held Dirac model with two major modifications: 1) there is a different 
dynamics for the baryons dictated by the lagrangian of eq. (H) and, 2 ) there 
is a clear way to uncover the baryon response to the medium through the 
soliton equation of motion. 

The many body skyrmion fluid in the dilute approximation is built from 
the product ansatz of eq. (|3|) for the skyrmions and the additive ansatz for 
the mesons of eq. (|4|). I 11 the following we use as input a distribution function 
for the baryons as a shortcut. 

A better method would be one that uses a fully self-consistent method, 
such as a Hartree-Fock. However, the use of a definite distribution function 





permits a first pass to the problem. We will use a Thomas-Fermi approxi¬ 
mation: The skyrmion fluid is forced to remain within the boundaries of the 
nucleus. This introduces a cutoff Fermi radius. Beyond this radius the meson 
fields will still exist, but the baryon densities will be taken as vanishing. 

We stress again that the results obtained here are to be taken within 
the context of the two approximations they hinge upon, namely, a simplistic 
product ansatz -dilute fluid approximation- and an averaging procedure with 
a pre-chosen distribution function. The former is a working hypothesis that 
permits the treatment of dynamics in a rather straightforward manner, and, 
the latter finds support in the the fact that we want to calculate equilibrium 
properties. The first approximation breaks down at short distances between 
nucleons and the second loses meaning for light nuclei, for which the ’mean- 
fields’ are not a wise choice of degrees of freedom, and direct nucleon-nucleon 
interactions are a better way to deal with the problem instead of by means 
of intermediaries. 

The skyrmions dynamics is obtained by boosting the static skyrmions 
rigidly. We perform a Lorentz boost on the collective coordinate R(t) of each 
skyrmion, although a non-relativistic Galilean transformation may suffice and 
perhaps be more consistent with the spirit of the rigid boosting. However, 
using a Lorentz transformation will help us to find the spin-orbit interaction 
below. 

For the sake of simplicity consider a Lorentz boost along the x axis with 
velocity parameter v (For a different view of the boost problem see ref||14j| ): 


x 


x 


y 


x — R(t ) 

y/1 — V 2 

y 


z = z 
F(r) -»■ F(r) 
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err 


o' r 


or 


B* 


uqrj vu{ r) 

a/1 — v 2 ’ a/1 — U 2 
^(r) vB 0 ( f) 


, 0,0 

, 0,0 


( 8 ) 


-\/l — u 2 ’ \/l — v 2 

with = uo(r)B 0 (r). 

Introducing the above transformation in eq. (]T|) and calculating the Hamil¬ 
tonian, we find the energy of a skyrmion in motion to be 

2 p 2 + 3 M 2 4 p 2 + 3 M 2 M 
f F 4 


E p — ( E2+E a — E iJ 
where 


3e M 


3e M 


^(U a -Uu+U in ^ (9) 


e = Jp 2 + M 2 , p = 


Mv 


\/l — v 2 


M is the skyrmion static mass for nonvanishing cuq, cto 


M(R) 
M(r ) 


r 00 

47t / r 2 drM ( r ) 

sin 2 F 


= e 


F 2 

2(7 ^ 7T 


and 


F' 2 + 2- 


+ 


1 sin 2 F 
2e 2 r 2 


sin 2 F 


+ 2F' 


i/2 


f 2 = 
f 4 = 

E a = 

Eu, = 

U a = 
U,„ = 




47 t f , ( »o sin 2 F 

2?h{ 2F + — 

—±fr 2 dr e 2 V 2 


47T J r 2 dr u/ 2 
47 t J r 2 dr V a 

47 r J r 2 dr e 2 "^ 
2gv 


2 , ,2 


[ dr u F' sin 2 F 
7T J 


( 10 ) 


( 11 ) 
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If the fields are kept as constant throughout the nucleon volume, we can 
use the skyrmion profile equation (or the virial theorem) to obtain 

M(R ) = e a M 0 (12) 

where Mo is the skyrmion mass for a = 0. We can include the rotational 
energy of the skyrmion in the adiabatic approximation. The result will yield 
the same scaling law as above. An attempt to include the rotational energy 
beyond the adiabatic approximation will induce the known instabilities to 
radiation of pions. Henceforth we use the measured free nucleon mass for M 0 , 
that contains both static and rotational energies. The above scaling is the 
only effect of the mean fields on the skyrmion, while the skyrmion contributes 
to the mean holds through its energy and baryon number separately. The 
topological character of the baryon current in the Skyrme model prevents it 
from changing. The only effect of the mean fields on the skyrmion is that of 
mass scaling. 

We can now write down the energy of N skyrmions in the mean held 
approximation with a Thomas-Fermi distribution 

/(Pi, R-i) = 6 (Pf(-R») ~ Pi) (13) 

with the normalization condition 

J dRi dpif = 1 (14) 

The Fermi momentum pp depends on the distance R, L of the center of the i th 
skyrmion from the center of the nucleus. Also, for indistinguishable particles 
Pf(Ri) = Pf(R)- The distribution function for the ensemble is just a product 
of functions of the type above. 

Averaging over the coordinates R and the corresponding momenta with 
the distribution function for zero temperature, the energy of the nucleus with 
spherical symmetry we obtain, using eq.(|J) 
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where 



E 


E(R ) 

E a 

= \l 

E w 

= —- u/ 2 - 

2 

I^int 

= g v u B 

Egk 

- 


J int 


E 


sk 


(15) 


,2 , ,2 


fk N 


M 2 


(16) 


where kp,k^ are the proton and neutron local (R dependent) Fermi mo¬ 
menta. 

In eq. fllGf) primes denote derivatives with respect to R } whereas in eq. (0) 
they represent derivatives with respect to r. Also 


B(R) = 


(kp) 3 (kprY 
3n 2 3 tt 2 


(17) 


The Euler-Lagrange equations for the mean holds become 


r 5 e 


2 a 


o '" + 2a' 2 + 


2cr' 

~R 


u" + 


dVa 2 2 2d- 9 E S k 

— + m u oj e - —— = 0 
da oa 

2u' 9 9 

— -mlui e 2 + g v B = 0 


(18) 


The ground state of the nucleus for fixed number of protons (P) and neu¬ 
trons (N) is obtained by minimizing the energy, constrained by Lagrange 
multipliers, with respect to the local wavenumbers kp and kpi 
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6E-fip 5P — n N 5N =0 
The algebraic equations for the multipliers become [|2j 


( 19 ) 


A*p — 9v w + 


g p b + q a + \Jkp + M 2 


I^N = 


v / br+ 9/ J 


( 20 ) 


In the Thomas-Fermi approximation, the nucleus has a finite radius be¬ 
yond which the nucleon densities are taken to vanish, therefore /ip = /xjv- 
In previous works JT^, [16|, we found that it is possible to fit the properties 
of nuclear matter, including the incompressibility, around 270 MeV, and the 
ground state densities of some magic nuclei, by introducing a parametrization 
of the scalar potential with 4 parameters. 

The key restriction for the terms in the extended potential comes from the 
requirement that the expectation value of the a field in the vacuum remains 
unchanged < uac|cr|uac >= 0. Terms of the form e na — 1 , generically referred 


to as no-log terms ||TT|, are acceptable. 

In the present work we will not commit ourselves as to the nature of the 
scalar field and its mass. We will not determine the scalar potential either. 
The field will be determined phenomenologically from the nuclear densities. 


3 Energy levels of a skyrmion in the nucleus 

In the present work we address the question of the energy levels and prop¬ 
erties of the nucleon in the skyrmion fluid. To this end we simplify the 
treatment of section 2. We solve the uj meson equation (|18|) with B(R) re¬ 
placed by the measured densities and the scalar determined by the chemical 
potential. In this manner we avoid tedious parametrizations of the scalar po¬ 
tential in order to fit the densities. The assumption is that such a potential 
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exists, and the results of ref. |L5| support it. The equation of the cr meson is 
therefore not solved. The uj meson eq. (|T^) is solved starting at the center 
of the nucleus up to a cutoff radius R c at which to=0. The a held is then 
extrapolated beyond that point using suitable boundary conditions^, [T(J. 

A key element of the nucleon interaction in the nucleus is the spin-orbit 
potential. We will here proceed to consider the contribution of both the 
scalar and vector fields to this force. The dominant contribution will arise 
from the modification of the mean fields in a rotating frame as viewed by the 
moving nucleon. The treatment will be a mixture of relativistic (when es¬ 
sential) and nonrelativistic methods. A more coherent treatment would call 
for a parametrization of the metric of a rotating nucleus, but this is outside 
the scope of the present work. 


Consider first a rigid boost of the skyrmion center with velocity v of eq. @ 
together with a rigid rotation of the skyrmion with collective coordinates 
decribing spin and isospin || 


U(i) -> A(t)U(r)A^ (t) (21) 

The spin-orbit interaction may then arise from the well-known Thomas 
precession, ft can be obtained from two consecutive Lorentz transforma¬ 
tions. This spin-orbit force will exist regardless of the extended nature of 
the skyrmion provided it moves in a central force held and it carries an axis. 
An easy way to implement this transformation is to take the isospin vector 
matrix as time dependent 


t = —f It x t 


( 22 ) 


where f 1 T is the Thomas frequency. Inserting the above ansatz in the skyrmion 
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lagrangian of eq. (|1) the spin-orbit interaction at the lowest order in the ve¬ 
locity is found to be 


U s . 0 . 


SL dV c 
~2 Ml R dR 


(23) 


where Vc is the central potential of the skyrmion in the nucleus, S is the 
spin and L the angular momentum. The same result as found in standard 
textbook derivations [[HJ]. In eq.(^3|) we have used the projection formula (3| 


i t A 


-ir-S 

2 A (!) 


(24) 


where A (R) is the moment of inertia of the nucleon[Q 


HR) 

A (R) 


[ r 2 dr A (R) 

o J 

sin 2 (F) 


2 „2cr 


Fi e 


4 /„, 0 sin 2 (F ) 


+ -\F' Z + 


(25) 


and F is the skyrmion prohle whose R dependence enters through the scalar 
field cr. It turns out, as expected, that the spin-orbit of eq. (^3|) is quite 
negligible, due to the -A? dependence. 

There is another source of spin-orbit interaction. It is due to the trans¬ 
formation of the fields to a rotating frame, essentially the coupling of the 
baryon current to the u held in a rotating nucleus analogous to the isoscalar 
coupling to the photon. However, the u meson coupling to the skyrmion, 
is not a gauge invariant one and we do not expect the same result as for 
the magnetic moment of the nucleon. The introduction of a gauge invariant 
coupling to the skyrmion requires additional terms. 

In order to find the spin-orbit potential, we need the mean fields for a 
streaming nucleus. In this case there arises a spatial component of the u 
meson held. An appropriate approximate ansatz for this component is []2| 


<jl> = V o;i(i?) = (fl x R) uq(.R) 


(26) 
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where V is the tangential velocity of the nucleus at each R and the an¬ 
gular velocity. At the same time the nucleon baryon density develops a time 
dependent piece of the form [[H| 


B(u) = S x u 


B 0 (u) 
2 A 


(27) 


where Bq(u) is the static baryon density of the skyrmion. 

Inserting eqs.(|26.27) above in the lagrangian there will appear new 


terms in energy. After averaging over the angular directions of R, we find 
the equation of motion for uq to be 


U " + ~R - + 3v B = 0 


(28) 


This is very similar to the equation of motion of the static to in eq. (h8|) . We 
solve equation (|28|) for each nucleus using the scalar held of the static case 
and demanding a vanishing uq at infinity. 

In order to find the corresponding spin-orbit interaction we consider a 
nucleon spinning at rest with a nucleus rotating with a velocity, — V(R), 
opposite to the direction of rotation of the nucleon. Using the collective 
coordinate quantization scheme of eq. O. and the projection formula of 
eq. (|24|) we find 

-SL 


Ws '°- 2 M 0 A (R) 


toi(R) 


(29) 


Clearly, W S . Q , is more important than U s , a , due to the -A- dependence. It 
is a pure skyrmion spin-orbit as evidenced by the presence of the moment of 
inertia in the potential. 

In the Dirac type of Walecka models |2|, the spin-orbit interaction arises 
from the coupling of the lower components of the Dirac wave function. In 
conventional nuclear forces calculations the spin-orbit interaction is found by 
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including pion exchange ||19|| . Here, it arises from the interaction of the rigid 
rotation of the nucleon with the flow of the mean fields, which is quite a 
different mechanism. 

The static fields contribute to the energy of the nucleus an amount pro¬ 
portional to the baryons rest energy, while the dynamical uq contribution 
originates from the rotational energy. The latter then adds a small correc¬ 
tion to the static mean fields and may be ignored. 

We now focus on the single skyrmion in the nucleus. Using eq.(||) in the 
relativistically invariant lagrangian of eq. the classical Hamiltonian of a 
single skyrmion in motion becomes 


H = V p 2 + M 2 + g v u; + W a . 0 . + U s . 0 . (30) 

where p is the nucleon momentum, the conjugated variable to the skyrmion 
center location R. Expanding the square root in eq. (|30|) to order p 2 and 
quantizing the coordinate R , we obtain an effective Schrodinger equation for 
the radial wave function of the skyrmion center with total energy E 


r d 2 2d 
.dR ? + RdR 


1(1 + 1 ) 

R 2 


Q(R) 


T = 0 


(31) 


where 


Q(R) = e 2 " Ml - ( - g v u + E - W,+ - Z(R ) 
Z(R) 


gv 


d 2 2d 

+ 


uj + 2 Ado U s . 0 . 


2M(R)ldR 2 ' RdR 
The central potential entering the calculation of U s . 0 . of eq.( 

Qi(R) 


V c = 


2 Mn 


(32) 
is given by 

(33) 
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where Q i is given by Q(R ) of eq.(p2|), but without the spin-orbit pieces. The 
energy eigenvalue of the Schrodinger equation E appears inside the operator 

Q{R). 

We have solved the Schrodinger equation ( fill) for the ground state single 
particle levels for the magic nuclei C 12 ,0 16 and Ca 40 . 

The only free parameters are the values of the scalar and the u meson 
at the center of the nucleus. These parameters determine the chemical po¬ 
tential and the cutoff radius. They are fixed once for all the energy levels in 
the chosen nucleus. We used the measured u meson and nucleon masses and 


g v = 7.31[1£ 


Table 1 shows the comparison between the 
the experimental ones averaged over proton 


predicted binding energies and 
and neutron states | 
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Table 1: Binding energies of single particle levels 


Nucleus 

Shell 

calculated energy 
MeV 

experimental energy 
MeV 

C 12 

lsi 

36.3 

35.2 


lp3 

2 

15.7 

16.9 

O 16 

lsi 

37.1 

43±5 


Ip | 

20.5 

20.1 


lpi 

2 

15.6 

13.9 

Co 40 

lsi 

48 

50±10 


lj>| 

35 

34±6 


lpi 

o 

30.7 

34±6 


Ids 

21 

18.5 


2sl 

15.7 

14.5 


1 dl 

2 

14.8 

12 
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The results show that the skyrmion picture of both the central and the 
spin-orbit interaction is quite good. The spin-orbit originates solely from 
the LV meson as viewed by the rotating skyrmion, in contradistinction to the 
Dirac mean field case in which both scalar and vector fields act together to 
produce a large interaction. It might appear at first that the skyrmion viewed 
as based upon a jj- expansion of QCD, cannot lead reliably to such fine detalis 
such as binding energies of around 20 MeV. However, the binding energies 
we find are not clearly related to the A- expansion. The A- approximation 
refers to the properties of the soliton constructed from colored quarks. Finite 
N c corrections, as well as quantum effects may indeed have influence on the 
basic properties of the soliton. However, the interaction of the skyrmion 
with the external scalar and vector fields operating through location of its 
center, promoted to the rank of quantum collective coordinate, is expected 
to be relatively unaffected; provided we use the appropriate parameters for 
the nucleon, namely its mass and moment of inertia. 

In the Skyrme model, the nucleon interacts as a colorless object with its 
neighbors. The ’exchanged’ meson fields are also colorless. Whether nucle¬ 
ons are partially deconfined inside the nucleus and quarks can percolate and 
be exchanged between nucleons, remains an open question. Nevertheless, in 
the language of skyrmeons, the leading contribution to the nucleon-nucleon 
interaction is indeed due to processes that ignore the color structure of the 
baryons, it is color-blind. This assumptions are based on the fact that the 
interactions between nucleons play a role in nuclear structure, only for in¬ 
ternucleon separations that are too large for color forces to be noticed. The 
effect of those forces is introduced phenomenologically in the parameters of 
the model, whose scale is not calculated. In some sense, the skyrmion in¬ 
teractions are like Van der Waals forces that act between neutral atoms and 
molecules. Although their origin resides within electromagnetic interactions, 
once the effective forces are determined one can work with them. The ap¬ 
proximation is even better for quarks because they are confined. 
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4 Nucleon properties in the medium 


In the previous sections we took the value of the meson fields as constant 


throughout the skyrmion. This method yields the scaling law of eq. (0). In 


this section we check the validity of this approximation and look for influences 
of the skyrmion interior on the meson fields and viceversa. Consider now the 
a and u fields radially symmetric, but allowed to change inside the skyrmion, 
namely 



w 0 = Wo (z) 


(34) 


where z is the distance from the center of the nucleus. Inside a skyrmion 
we have 



(35) 


with R the center of mass location of the skyrmion and r the radial 
distance from this center. 

We will simplify the calculation by averaging the fields over angles. The 
distribution function we use is independent of angle (eq. ( |I3|) ). The soli- 
ton profile will still be taken as spherically symmetric. We could allow for 
distortions of the nucleon, but this will make the calculation extremely cum¬ 
bersome. We want to keep the description as simple as possible and still 
capture the esence of the effects of the meson fields on the skyrmion. 

We average over angles, without determining the meson fields yet 
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P(r, = (4^)2 / / w(z) 


(36) 


With x = e 2fJ obtaining 


| rR+r 

Q( r i R ) = ~ A —f/ x(u)udu (37) 

4 rP 7|H-r| 

with a similar equation for P. 

The single skyrmion static energy is now 


/* sin 2 F F' 

E = 4n J r 2 dr (Q(r, R) M(r ) - P(r, P)gy ^ 2 ) (38) 

where M(r) is defined in eq. ([10]). The meson fields depend in a nontrivial 
way on R and the skyrmion inherits this dependence. 

The skyrmion equation of motion becomes 


/ 8 sin 2 F 


F" 


4 sin ‘IFF' 1 Q sin 2F 


+ F'(2 Q/r + Q') + _ 2 

4sin 2 F sin 2F 2 g v P' sin 2 F 

= 0 




7r 2 r 2 


(39) 


where r = eF n r and primes denote derivatives with respect to r. 

We start the calculation by evaluating the meson fields for the chosen 
nuclei. The values of the meson fields at the center of the nucleus are varied in 
order to fit the nuclear densities. We do not parametrize the scalar potential 
in a definite manner. We take the measured densities as input and solve the 
Thomas-Fermi equations of motion for the u meson with a fixed chemical 
potential /j. The meson fields are then fitted with analytical expressions that 
reproduce accurately the calculated values. For the nuclei C 12 and O 16 we 
use formulae that are easily integrable in order to find Q and P of eq. (|3G|) 


22 










Figure 1: M e ff in C 12 , constant (modified) meson fields inside the skyrmion, 
full line (dashed line) 


x(z) = 1 + (c?i + bi z 2 + ci z 4 ) e dl z2 

u{z) = (a 2 + b 2 z 2 + c 2 z 4 ) e~ d2 (40) 


whereas for, Ca 40 , it was found after trial and error that a better suited 
expression for the u held is 


u(z) = (a 2 + b 2 z 2 + c 2 £ 4 ) e“ d2 2 


(41) 






R (fm) 


Figure 2: M e ff in O 16 , constant (modified) meson fields inside the skyrmion, 
full line (dashed line) 


We then calculate the new soliton profiles using eq.([39|), evaluate the 
skyrmion effective mass using eq. (|38|) and compare it to the desired effective 
mass needed in order to fit the single particle levels obtained in the previous 
section as good as possible. We iterate the procedure until the effective 
masses are the closest to the desired values. 

Table 2 shows the parameters for the analytical expressions of eqs. 


W- 


Figures 1 through 3 compare the obtained effective skyrmion masses both 
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R (fm) 

Figure 3: M e ff in Ca 40 , constant (modified) meson fields inside the 

skyrmion, full line (dashed line) 

with and without meson fields variation inside the skyrmion. 

The effective mass of the skyrmion is defined as Mo e a ° for the case of a 
constant scalar field inside the nucleon, and 

Meff = E + T (42) 

where E is the nucleon energy obtained from eq. ( |38| ) and A is defined in 
eq. ( |25|) . The trend and absolute values of both types of effective masses are 
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Table 2: Parameters for the fits to the meson helds with the analytical 
expressions of eqs. [42,43] 


parameter 

C 12 

O w 

Ca 4U 

Gp 

-0.4744275 

-0.3581017 

-0.5112333 

bi 

-0.2732692 

-0.2322270 

-0.1049903 

Cl 

-7.1264720E-03 

-2.4487255E-03 

-1.7629359E-02 

d\ 

0.4292937 

0.3366868 

0.2366882 

a 2 

0.1471073 

0.1004175 

0.1679764 

b 2 

5.3953879E-02 

5.9962720E-02 

0.2383440 

c 2 

0.1121303 

2.7225710E-02 

-9.0451920E-03 

d 2 

0.8262799 

0.5333008 

0.9995874 


consistent. However, they are not identical. It was impossible to improve the 
agreement beyond the one shown in the graphs. The inclusion of the volume 
of the skyrmion in the evaluation of the meson helds cannot be disregarded. 
This is more clearly seen in figure 4, where we show the scalar held x — c 2a 
in the nucleus with and without the influence of the volume of the skyrmion 
for O 16 . A similar picture emerges for the uj held and in the other nuclei. 

Another dramatic effect shows up in the isoscalar root mean square radius 
of the nucleon. Figure 5 depicts this radius as a function of the location of 
the center of the skyrmion in O 16 . The nucleon swells inside the nucleus to 
almost twice its free size. This is a bit troubling, because it would imply a 
breakdown of the fluid approximation. The nucleons no longer act as free 
molecules and the solitons overlap. However, the overall picture does not 
change because the strength of the interaction between solitons dimishes 
accordingly. A measure of that interaction is the 7 tNN coupling constant 0]. 
A better estimate of the importance of the overlap between skyrmions will 
demand a calculation of the skyrmion-skyrmion interaction, but we limit 
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R(frn) 


Figure 4: y as a function of R in O 16 for constant (modified) meson fields 
inside the skyrmion, full line (dashed line) 

ourselves to obtain a rough estimate of the effect. 

Figure 6 shows this constant as a function of distance for O 16 . It has 
the very opposite behavior. It drops to almost half its value inside the nu¬ 
cleus. So, perhaps there is still room to consider the nucleons as essentially 
free objects. Despite the fact that they overlap, the strength of the mutual 
interaction renders this overlap quite ineffective. 

Other nucleon observables, such as magnetic g-factors, axial coupling con¬ 
stant and magnetic radii are less susceptible to changes in the environment. 
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Figure 5: Skyrmion rms radius in O 16 with modified meson fields 

We have attempted to put the scalar fields we found in correspondence 
with the Kisslinger type of potential with the identification ||22||, 

—k 

X = ~ , 

1 + g 0 n 

n = An cop(R) (43) 

where p is the nuclear density, cq ~ .21 m~ 3 , and g' 0 < 1. This procedure 
would predict the values of these parameters as suggested by the Skyrme 
model. Although it was possible to fit the scalar field with the above potential 
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Figure 6: g n NN in O 16 with modified meson fields 

around the surface of the nucleus, it is impossible to obtain a moderately 
reasonable fit inside the nucleus. In all cases there seems to be convergence 
to a set of parameters that indicate a much larger Landau parameter than 
expected, approximately around g' {) ~ 1. The Skyrme model appears to be 
suggesting a strong medium correction to the so-called effective diffraction 
index. Inside the nucleus nonlinearities play a dominant role. 
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